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Abstract 

A thermodynamic model describing phase transitions with thermal memory, 
in terms of an entropy equation and a momentum balance for the microforces, is 
adressed. Converge results and error estimates are proved for the related integro- 
differential system of PDE as the sequence of memory kernels converges to a multiple 
of a Dirac delta, in a suitable sense. 
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1 Introduction 

In this paper, we deal with the following integrodifferential PDE system, describing a 
phase transition process in the case when thermal memory effects are included. Indeed, 
here d T stands for the absolute temperature, X T for the phase parameter, and k T for a 
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memory kernel 



<9t(lntf T + X(X T )) - k q M t - A(k T * tf r ) = / (1-1) 

d t Xr-AX r + Z T + <T , (Xr) = \'(X T )& T , £ r G /3(X T ) (1.2) 

$ r | r = $ r and <9„X r | r = (1.3) 

lntf r (0) = ln$ and X T (0) = X . (1.4) 



Each of the partial differential equations fll.ip — (jl.2p is meant to hold in a three-dimensional 
bounded domain Q, endowed with a smooth boundary T, and in some time interval (0, T). 
In f ll.ip . the memory kernel k T may depend on a positive parameter r. Moreover, the sym- 
bol * denotes the usual time convolution formally defined by (a * b) (t) = J * a(t — s) b(s) ds 
for functions that depend just on time, and then extended to functions that also depend 
on space. Furthermore, / is some given source term. In (11.21) . (3 is a maximal monotone 
graph in M 2 , while A and a are real functions defined on the whole of 1R. The bound- 
ary conditions (jl.3p must be satisfied in T x (0,T), while the initial conditions ( II. 4p are 
written for the functions lni? T and X T : of course, $o, and Xo are given boundary and 
initial data. 

Equation fll.ip may be interpreted as an entropy balance equation. Note in particular 
that that the equation is singular with respect to the temperature, mainly for the presence 
of the logarithm, forcing the temperature to assume only positive values (which is in 
accordance with physical consistency). Similar systems have been studied in the literature 
from the point of view of the existence and regularity of solutions (see, among the others, 

[si \mu u cms]). 

The well-posedness of a proper variational formulation of ( ll.ip - (jl.4p has been proved 
in [3] . Here, our main goal is the following. By assuming that k T converges to n' 5 at r \ 
in a suitable sense, where 8 is the Dirac mass at the origine of the real line and k' is a real 
constant satisfying k := k + k' > 0, we prove that the solution (i? T ,X r ) to (jl.lD — (jl.4p 
converges in a proper topology to the solution ($, X) of the problem stated below 



d t (\n$ + \(X)) -«A0 = / (1.5) 

d t X- AX + £ + a'(X) = \'(X)d, £eP(X) (1.6) 

■&\ T = ^ r and d„X\ r = (1.7) 

lntf(0)=ln$ and X(0) = X . (1.8) 



This convergence result is obtained by the use of an a priori estimates technique and 
passage to the limit arguments, based on monotonicity and compactness. Moreover, an 
error estimate, i.e. an estimate of suitable norms or quantities involving the difference of 
solutions, is shown. 

Our paper is organized as follows. In the next section, we discuss a derivation of 
the system (ll.5p - (ll.6l) from the basic laws of thermomechanics. Section [3] is devoted 
to the statement of our assumptions and of our results on the mathematical problem. 
In Section HI we present some auxiliary material that is needed for the proof of our 
convergence Theorem I3.3| mainly. The last section is devoted to the proofs of the above 
theorem and of the error estimate stated in Theorem 13.41 
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2 The model 

In this section, we briefly introduce the modeling derivation of the equations fll.ll) -f TL~2]) 
and discuss the convergence to ( 1 1 . 5f) — ( TTT6T) . as the parameter r (in the memory kernel) 
tends to 0. Here, the argument is mainly developed from a physical point of view, while 
we refer to subsequent sections for a more precise setting of analytical assumptions and 
comments. In particular, we aim to focus on the fact that (11.11) accounts for thermal 
evolution involving memory effects, on the basis of the memory kernel k T . 

Materials with thermal memory have been deeply studied in the literature, both from a 
modeling and analytical point of view. We refer, in particular, to the approach by Gurtin 
and Pipkin (see [15]) for thermal memory materials. Several authors have investigated 
phase transitions in special materials with thermal memory, both concerning modeling 
and analysis. For a fairly complete and detailed presentation of this kind of problems, 
let us mention the very recent monograph [jQ. Now, we combine thermal memory with 
a new theory for phase transitions models, based on a generalization of the principle of 
virtual powers (see [H]). The idea is that micro-forces, which are responsible for the phase 
transition, have to be included in the whole energy balance of the system. Consequently, 
the phase (evolution) equation is derived as a micro-forces balance equation and it is 
coupled with an entropy evolution equation. This approach has been recently investigated 
in the literature by several authors (among the others, we mainly refer to the papers [3] 
and [7J, in which the derivation of the model is detailed in the case when possible thermal 
memory effects are included, as in equations ( ll.ip -( TL~2]) ). 

Indeed, let us recall that in [3] the theory by Gurtin-Pipkin is considered, allowing 
the free energy functional to depend on the past history of the temperature gradient. 
The resulting functional accounts for non-dissipative contributions in the heat flux, which 
may be combined with additional dissipative instantaneous contributions coming from a 
pseudo-potential of dissipation. The use of an entropy balance has been recovered, in 
this approach, from a rescaling (with respect to the absolute temperature) of the energy 
balance, under the small perturbations assumption (see also |H [5]). In [7] a fairly general 
theory is introduced. The model is derived by a dual approach (mainly in the sense of 
convex analysis) in which the entropy and the history of the entropy flux are chosen as 
state variables (together with the phase parameter and possibly its gradient). Then, the 
dissipative functional is written in terms of a dissipative contribution in the entropy flux 
and for the time derivative of the phase parameter. 

Let us point out that the above mentioned approach is not far from the theory proposed 
by Green-Naghdi [13] and Podio-Guidugli [16], in which some thermal displacement is 
introduced as state variable (it is a primitive of the temperature) and the equations 
come from a generalization of the principle of virtual powers, in which thermal forces are 
included. As a consequence, in this framework, the entropy equation is formally obtained 
as a momentum balance (i.e., a balance of thermal forces acting in the system). The 
reader may also examine [HI [10] , where some asymptotic analyses are carried out to find 
the interconnections among different Green and Naghdi types. 

We aim to observe that the model we are investigating actually may be obtained 
by combining the above two theories, i.e. generalizing the principle of virtual powers 
accounting for microforces as well as thermal stresses. Let us present our position. First, 
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we specify the expression of the power of internal forces. The power of interior forces is 
written for any virtual micro- velocity 7 and thermal velocity v, as follows 

3>i= / 5 7 + H- V 7 + Q- Vv, (2.1) 

where B and H are interior forces responsible for the phase transition (as introduced 
in [H]), and Q stands for a thermal stress (corresponding to the entropy flux by [16]). 
Hence, the resulting balance equations are written as momentum balance equations. It is 
assumed that an external (density of) entropy source / is applied. A thermal momentum 
is introduced to measure reluctance to the order of the system (in analogy with the 
mechanical momentum measuring reluctance to quiet). We prescribe that it is given by 
the entropy s. It results that (see (12.11) ) 

s t + div Q = /. (2.2) 

As far as the microscopic momentum balance, we assume that no acceleration and no 
external forces are contributing, so that we have 

B — div H = 0. (2.3) 

Henceforth, fl2.2l) - fl2.3p are combined with suitable boundary conditions. As usual, we 
assume that the flux through the boundary H ■ n is null, while (mainly for analytical 
reasons) we prescribe a known temperature on the boundary. 

The entropy s, the entropy flux Q, and the new interior forces B and H are recovered 
by suitable energy and dissipation functionals, that we are going to make precise, in terms 
of state variables. The state variables are related to the equilibrium of the thermodynam- 

ical system: they are the absolute temperature the phase parameter X, the gradient 

• — - 1 

VX (actually accounting for local interactions), and the history variable V$ , which is 
defined as 

W (s) = / W(r)dr, s > 0. (2.4) 

Jt-S 

I— I — * 

As in [3], we assume that the free energy of the system (depending on V$ , X, VX)) is 

split into two contributions: the first is related to present variables at time t (^p), the 

second accounts for some history in the system (^h), measured through a memory kernel 

(related to k T in the equations). In particular, the history contribution of the free energy 

is given by 

= i|W* H T (2.5) 

where S T is the space of the past histories (as it is introduced in the theory of thermal 
memory materials by Gurtin and Pipkin), defined by 

r+00 

S r := {f : (0, +00) -> M 3 measurable s.t. / h T (s)\f {s\ 2 ds < +00}. (2.6) 

Jo 

Here, h T : (0, +00) — > (0, +00) (possibly depending on a parameter r) is a continuous, 
decreasing function such that 



r+00 

I s 2 h T (s)ds < +00. 




(2.7) 
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The space S T is endowed with the natural norm 



|2 



+00 



h T (s)\i(s)\ 2 ds (2. 



and the related scalar product is (v, u) St = J +o ° /i r (s)v(s) ■ u(s)ds. Let us comment that 
in our system, to derive (11.11) . we have introduced a kernel k T such that —k' T = h T . More 
precisely, let k T : (0, +00) —> R and require that 

fere^^+oo), \imk T (s) = 0. (2.9) 

s— >+oo 

Hence, by virtue of the assumptions on h T we also have 

k' T < and > a.e. in (0, +00). (2.10) 

Note that k' T (t) vanishes for t going to +00 and that k T is a non- increasing function with 
^r(O) > 0, and in the case k T (0) = one has k T = 0. These assumptions on k T actually 
ensure that the model is thermodynamically consistent, as it is detailed in [3J. 

Then, the free energy functional typ (written at the present time t) is addressed 

* P (i?, X, VX) = c v $(l - lntf) - A(X)0 + <r(X) + /3(X) + i|VX| 2 (2.11) 

where cy > (in the sequel let us take cy — 1) is the specific heat, a and A are suffi- 
ciently smooth functions (with A'(X) denoting the latent heat), (3 is a proper convex and 
lower semicontinuous function, possibly accounting for internal constraints on the phase 
variable X. For instance, a fairly classical choice is /3(X) = 7[o.i](X), which is equal to if 
X G [0, 1] and takes value +00 elsewhere (thus forcing X £ [0, 1]). 

Dissipation is rendered in terms of the time derivative X t and of the dissipative vari- 
able V$. It is derived by a pseudo-potential of dissipation (in the sense by Moreau, i.e. 
a convex, non-negative function assuming its minimum for null dissipation): 

$(X f ,W) = ^|X t | 2 + ^|W| 2 . (2.12) 

Note that, in order to ensure the validity of the second principle of thermodynamics, it is 
required that k > 0. 

Now, we are in the position of recovering our system, after specifying constitutive 
relations for the involved physical quantities. We have that 

s = -|| = lntf + A(X) (2.13) 

and 

B G w + H = m + a ' {x) ~ + Xt {2M) 

(3 being the sub differential (in the sense of convex analysis) of /3, and 



(2.15) 
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Hence, the entropy flux vector Q is specified by 

- Q = {-Q nd ) + {-Q d ) (2.16) 

where — Q nd results to be defined in § T . It is obtained taking the derivative in S T of the 
history functional with respect to the history variable. Integrating by parts in time, using 
the Frechet derivative, and exploiting the hypotheses on k T (see [5] for any further detail) 
lead to 

- Q nd (t) = - f k T (t - s)W(s)ds. (2.17) 
We have now to make precise the dissipative part of the entropy flux 

-Q d = ^ = KqW. (2.18) 

Hence, equations fll.ip - fll.2p are obtained by (12. 2 p and (12. 3 p exploiting the above 
introduced constitutive relations. We point out that in (II. ip . the past history contribution 
of (12.171) (actually its divergence), i.e. J_ k T (t — s)'Vi)(s)ds, is assumed to be known and 
included in the external entropy source / (we have used the same notation of (12. 2p for 
the sake of simplicity). 

As we have already pointed out in the Introduction, the main aim of this paper is to 
investigate the asymptotic behavior of system ( Il.ip -( 1L~2"|) as the thermal memory kernel 
converges to k' 5, 5 being the Dirac mass at the origin of the real line and k' > 0, in a 
suitable sense 

k T k' 5. (2.19) 

More precisely, we are interested in proving that solutions to the system (ll.ip -( TL~2]) 
converge to solutions to f ll.5p -f fL6|) (at least in some weak topology). Let us briefly 
comment that the system fll.5p - fll.6p . obtained in our proof as a suitable limit of (11. ip - 
(11. 2p . can be actually derived by an analogous procedure as the one we have performed to 
formally derive (JHI])-(jL2l). Indeed, (TL5l)-(Tr6l) follow from $2J|) and (E| when exploiting 
(I2.13p -( l2.16p . Here, the new energy and dissipative functionals are \1/ = \l/p (i.e., no history 
contribution of type (12. 5p in the free energy is given) and (cf. (I2.12p ) 

9(Xt, W) = + (^ + ^o) |w|2- (2 _ 20) 

In particular, it results that in fl2TT6|) Q nd = 0, while (due to (f2TT8"]l ) -Q d = (k + «' ) W. 



3 Statement of the mathematical problem 

In this section, we make our assumptions precise and state our results. First of all, we 
assume Q to be a bounded connected open set in M 3 (lower-dimensional cases could be 
considered with minor changes) whose boundary V is supposed to be smooth. Next, we 
fix a final time T G (0, +oo) and set: 

Q := n x (0, T), S := T x (0, T) (3.1) 
V'.= H\Q), V Q :=Hl{Sl), H:=L 2 (n) (3.2) 
W := {v E H 2 (n) : d„v = on T}, (3.3) 



Bonetti Colli Gilardi 



7 



d v denoting the normal derivative operator on the boundary. We endow the spaces (j3.2|) - 
(13. 3p with their standard norms, for which we use a self-explanatory notation like || • \y. 
Moreover, for p G [1, +oo], we write || • || p for the usual norm in L P (Q); as no confusion 
can arise, the symbol || • || p is used for the norm in L P (Q) as well. In the sequel, the same 
symbols are used for powers of the above spaces and the corresponding natural induced 
norms. It is understood that H C Vq* as usual, i.e., any element u G H is identified with 
the functional Vq 3 v \-> f n uv which actually belongs to the dual space V^* = iJ~ 1 (fi) 
of Vq. We observe that L 2 (0,T; V^,*) coincides with the dual space of L 2 (0,T; Vo) and use 
the symbol ( • , • ) for the corresponding duality pairing. 

As far as the structure of the system is concerned (see (11.11) . (11.51) and (11.21) . (11.61) ). we 
are given the three functions (3, A, a, the constant hq and the memory kernel k T depending 
on the parameter r > and we assume that the conditions listed below are satisfied. 



(3 : 1R — > [0, +oo] is convex, proper, lower semicontinuous, and /3(0) = (3.4) 

A, <j G C 1 (M) and A', a' are Lipschitz continuous (3-5) 

k >0 and k T eW hl (0,T). (3.6) 

We define the graph (3 in M. x R by 

13 := 3/3 (3.7) 



and note that (3 is maximal monotone and that /3(0) 3 0. In the sequel, we write D(f3) 
and D((3) for the effective domains of (3 and (3, respectively, and we use the same symbol 
(3 for the maximal monotone operators induced on L 2 spaces. 

As far as the data of our problem are concerned, we assume that the functions /, 



■$0; an d the constants and 0* are given such that 

feL\0,T-H) (3.8) 

tfr G H l (0, T; H" 1 / 2 ^)) n W^(0, T; L°°(r)) n C°([0, T}-H l ' 2 (T)) (3.9) 

< 0» < <&* < +oo (3.10) 

0* < #r < a.e. on Y x (0, T) (3.11) 

tfo e L°°(Q), $*<$ < a.e. in (3.12) 

X G V and /3(X ) G L 1 ^). (3.13) 



The function $r is the boundary datum for the temperature and we would like to consider 
a function u := d—d^x vanishing on the boundary as associated unknown function. Hence, 
a natural choice of is the harmonic extension of so that Am = A-#. Therefore, we 
define : Q — > R as follows 

0«(t)eV, A0 M (t) = O, and Mt)\r = Mt) for a.a. t G (0, T). (3.14) 

For the regularity of t?^ (induced by (13.91) ) see the subsequent Proposition 14.11 

Next, we list the a priori regularity conditions we require for any solution ($, X, £) of 
either (Q1)-(IL"4T) or (H^-JTB}. We ask that 

tfGL 2 (0,T;V0 and u := - ^ G L 2 (0, T; Vb) (3.15) 

0>O a.e. in Q and In G if 1 (0, T; V^") D L°° (0, T; H) (3.16) 

X G H^O, T; if) n L°°(0, T; 1/) n L 2 (0, T; W) (3.17) 

£gL 2 (0,T;#). (3.18) 
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At this point, we are ready to state the problems we are dealing with in a precise form. 
For fixed r > 0, we look for a triplet (i? T , X T , £ T ) satisfying fl3.15p - fl3.18p and the following 
system 

d t (lnA(t) + A(X T (t))) - K AMt) ~ A(fcr * #r)(*) = /(*) 

in Vq*, for a. a. t G (0,T) (3.19) 
d t X T - AX T + £ r + a'(X T ) = X'(X T ) A and £ T G 0(X T ) a.e. in Q (3.20) 
(lntf T )(0) = ln^o and X r (0) = X . (3.21) 

We note that the boundary conditions fll.3p are contained in (13.151) and (I3.17P (see the 
definitions f)3.14p and fl3.2p - fl3.3p ). We also remark that (I3.16P implies that ln-$ is a 
continuous Vq 1 - valued function (while no continuity of $ is known), so that the Cauchy 
condition for ln$ T contained in f)3.2ip makes sense. Similar remarks hold for the limit 
problem we are going to state (i.e., fll.5l) - fll.8p in a precise form). 

The following well-posedness result deals with a fixed r > and essentially follows 
from [5]. Just the notation is different, indeed. 

Theorem 3.1. Assume that both (13.41) -f l3~7l) and f l3.8p -f l3~T3j) hold. Then, there exists a 
unique triplet X r , £ r ) satisfying fl3.15p - fl3.18p and solving problem fl3.19l) - fl3.2ip . 

Our aim is to study the limit of the solution ($ T ,X T ,!; T ) as r tends to zero, under 
suitable assumptions on the behavior of the memory kernel k T . Namely, we assume that 

1 * k T k' strongly in L\0, T) as r \ (3.22) 



for some real constant k' q and set 



K + 4 . (3.23) 



In (I3.22p and later on, we use the same symbol for any real constant (like 1 and n' ) and 
for the corresponding constant function. We advice the reader that k > in the sequel, 
either by assumption or as a consequence of some condition we require, so that the limit 
problem we are going to state is parabolic with respect to Such a problem consists in 
looking for a triplet (■#, X, £) satisfying fl3.15l) - fl3.18l) and the following system 

d t (\n$(t) + \(X(t))) -«A0(i) = fit) in V *, for a.a. t e (0,T) (3.24) 
d t X - AX + f + a'{X) = \'{X) $ and £ G (3{X) a.e. in Q (3.25) 
(lni?)(0) = ln^o and X(0) = X . (3.26) 



By just taking k T = and replacing kq by k in Theorem 13.11 we obtain 

Corollary 3.2. Assume that both fETD -fl3T7l) and fl3~8l) - fl3~T3|) hold and that k>0. Then, 
there exists a unique triplet (-d,X,^) satisfying f)3.15l) - fl3.18p and solving problem f)3.24p - 
(J32SD. 

However, we can prove a convergence result under further assumptions, namely 
/ (K v(t) + (k T * v)(t))v(t) dt > 11^11^2(0^) and || A; r || L i (0 T) < n* (3.27) 
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for some constants k*,k* > and every v G L 2 (0,T) and r > 0. By taking v = 
on (T',T), we clearly see that the time T can be replaced by any T' G (0,T) in the first 
inequality of (I3.27p . Moreover, we observe that f)3.22p and (I3.27P imply that n > k*, so 
that k > as a consequence. Here is our first result. 

Theorem 3.3. Assume f l3.4p -f l3~T|) . (I3.8p - fl3.13p . and let ($ r , X T , £ T ) be the unique solu- 
tion to problem fl3.19p - fl3.2ip given by Theorem \3.1\ Moreover, assume f l3.22p -f l3.23p 
and (I3.27p . Then, ($ r ,X r ,£ T ) converges in a proper topology to the unique solution 
to problem fl3T2^) -f lH^j) satisfying f l3TT5|) - fl3~TH|) . 

The topology mentioned in Theorem 13.31 will be clear from the proof we give in Sec- 
tion 13.31 and is rather strong. Provided that a much weaker topology is considered, an 
error estimate can be proved. We have indeed 

Theorem 3.4. Under the assumptions of Theorem VS.'A the following estimate holds true 

||1 * (l? r — $)||ioo( 0iT;V /) + \\X T — ^||| O o( 0>T; ^) nL 2( 0iT; y) 

+ / (lntf T -ln0)(0 T -0)+ / (£t-0(*t-X) 

J Q J Q 

<M\\l*k T -K' \\ LHOtT) (3.28) 

for r small enough, where M depends on the structure and the data, only. 

Remark 3.5. We observe that a sufficient condition for 03.271) 9 is that k T has the form 
k T (t) = T~ l k{t/r) with k G -^ 1 (0, +oo). Moreover, a sufficient condition for (13. 271b is that 
k T is a positive type kernel, i.e., 

T 

(k T *v)(t)v(t) dt > for every v G L 2 (0,T). (3.29) 

In such a case, one can take «;* = Kq, indeed. The fact that the kernel is of positive type 
is actually in accordance with the assumptions required on k to ensure thermodynamical 
consistency of the model [3]. We remark that a sufficient condition for a kernel k to be 
of positive type is the following (see, e.g., [21 Prop. 4.1, p. 237] or [H]): k is smooth, 
nonnegative, decreasing, and convex. So, any positive multiple of exp(— t/r) plays the 
role and the kernel given by k T {t) := (k' /t) exp(— t/r) is a prototype for both f)3.22p 
and f!3.27p since 1 * k T converges to k' strongly in L 1 (0,T). More generally, assume that 
k T can be split as k T = p T + r T , where p T is of positive type and r T is a remainder. If 

1 * p T — y k\ and 1 * r r — > k! 2 strongly in L 1 (0, T) 
|KHli(o,t) < Vo for every r > 

for some constants k[, k' 2 , and r] < k , then both assumptions f)3.22p and (I3.27j) are still 
fulfilled. Indeed, we can take clearly, and k, = k — rj , since 

/ (kqv + k T * v)v dt > (k v + r T * v)v dt > (k — rj ) / v 2 dt for every v G £ 2 (0, T) 
Jo Jo Jo 

thanks to the Holder and Young inequalities (see also (I3.37P ). Note that k[ > since p T 
is of positive type, while k' 2 can be any real constant. 
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Remark 3.6. Assumption ( I3.22p is a well-defined reinforcement of the condition roughly 
mentioned in the Introduction as k T — > k' 5, where 5 is the Dirac mass at the origin. 
Indeed, if we introduce the Heaviside function H on (— oo, T), i.e., H(t) = for t < and 
H(t) = 1 for t £ (0,T), and the trivial extension k T of k T , (I3.22p reeds 

H * k T — >■ hi' H strongly in oo, T) 

with an obvious new meaning of the convolution. By differentiating and observing that 
(H * k T )' = 5 * k T = k T , we deduce that 

k T — > k' 5 in the sense of distributions on (— oo,T) 

where 5 is the actually well-defined Dirac mass at in the open set (— oo, T). 

Remark 3.7. By checking the proofs in the next sections, the reader will be able to 
realize that our results can be suitably extended to the case of coefficients kq t possibly 
depending on r, with boundedness and convergence properties as r \ 0. 

We recall that Q is bounded and smooth. So, throughout the paper, we owe to some 
well-known embeddings of Sobolev type, namely V C L P (Q) for p £ [1,6], together with 
the related Sobolev inequality 

IMIp — f° r every v £ V and 1 < p < 6 (3.30) 

and W l ' p (VL) C C°{Tl) for p > 3, together with 

IMloo < CplMlw^n) f° r every v £ W l,p (Q) and p > 3. (3.31) 

In (I3.30p . C depends only on Q, while C p in (I3.3ip depends also on p. In particular, the 
continuous embedding W C W 1,6 (Q) C C°(Q) holds. Some of the previous embeddings 
are in fact compact. This is the case for V C L 4 (f2) and W C C°(Q). We note that also 
the embeddings W C V, V C H, Vq C H, and H C Vq are compact. Moreover, we often 
account for the well-known Poincare inequalities 

\\v\\ v < C\\Vv\\ H for every v £ V (3.32) 
\\ v \\v < c(\\ Vf \\h + Un^l) for every f £ ^ (3.33) 

where C depends only on Q. Furthermore, we repeatedly make use of the notation 

Qt := Q x (0, t) for t £ [0, T] (3.34) 

and of well-known inequalities, namely, the Holder inequality and the elementary Young 
inequality: 

ab < 5a 2 + — b 2 for every a, b > and 5 > 0. (3.35) 

As far as properties of the convolution are concerned, we take advantage of the elementary 
formulas (which hold whenever they make sense) 



a *b = a(0)(l * b) + a t * 1 * b and (a * b) t = a(0)6 + a t * b 



(3.36) 
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and of the well-known Young theorem 



U * v\\ L r( 0tT . X ) < \\u\\lp(0,T)\\v\\li(0,T;X) 



(3.37) 



where X is a Banach space, 1 < p, q, r < oo, and 1/r = (1/p) + (l/<?) — 1 (cf., e.g., [H]). 
Finally, again throughout the paper, we use a small-case italic c for different constants, 
that may only depend on Q, the final time T, the shape of the nonlinearities A, /3, a, and 
the properties of the data involved in the statements at hand; a notation like c$ signals 
a constant that depends also on the parameter 5. The reader should keep in mind that 
the meaning of c and c$ might change from line to line and even in the same chain of 
inequalities, whereas those constants we need to refer to are always denoted by capital 
letters, just like C in (I3.30p . 

4 Auxiliary material 

This section contains a very short summary on the properties of the harmonic extension 
d'x of the boundary datum -#r (see (13.141) ) and a preliminary result dealing with a gen- 
eralized version of the limit problem fl3.24p - fl3.26p . The properties listed in the following 
proposition will be extensively used in the sequel. 

Proposition 4.1. Assumptions fl3.9p - fl3.lll) yield 



More precisely, owing to the theory of harmonic functions, in particular to the maximum 
principle, we have that 



where C is a constant depending on f2 ; only. 

Now, in order to help the reader, we sketch the outline of the proof of Theorem 13.31 
we are going to develope in the next section. By accounting for a number of a priori 
estimates and using well-known compactness results, we derive that the family of solutions 
(■# r; X r , £ r ) converges (for a subsequence) to a generalized solution to problem (13.241) - 
f)3.26p . in which ln-$ is understood in a non standard sense. Next, in order to conclude 
that such a solution actually is the solution given by Corollary I3.2[ we prove a preliminary 
well-posedness result for generalized solutions (Theorem 14.21) . Therefore, we first have to 
introduce the ingredients that are needed to explain such a notion of solution. 

We define a generalized logarithm by following [12], Def. 4.2]. However, we confine 
ourselves to consider the functions $ G L 2 (0, T; V) such that = t?r on the boundary, 
i.e., $ G '&-K + L 2 (0,T; Vq). First, we introduce the function -0 : M — > (— oo, +oo] by setting 



G H 1 (0,T;H)DW 1 ' 1 (0,T;L co (n))r]C ([0,T};V). 



(4.1) 



||^M|U 2 (o,T ; y) < C||^r||L2( 0T;H i/2( r )), 

II^mIU^O.T;^) < C||'^r||Li(0,T; J H'- 1 /2(r)) ; 

<C||^ r || 
< $x < a.e. in Q, 

||^^M|U 1 (0,T;L°°(n)) = 1 1 ^t^r 1 1 (0,T;i°° (r) ) 



ip(r) = r(logr - 1) if r > 0, ip(0) = 0, and i[)(r) = +oo if r < 0. (4.2) 



12 Singular limit of a integrodifferential system 



Then, for •& G $x + L 2 (0, T; V ), we term Ln$ the set of ( G L 2 (0, T; V*) satisfying 

(C, v - ■&} + [ < / ip(v) for every v G d + L 2 (0, T; Vb) (4.3) 

where ( • , • ) stands for the duality pairing between L 2 (0, T; Vb*) and Z/ 2 (0, T; Vb). It can be 
checked (see [121 Thm. 4.7])) that d > a.e. in Q whenever Ln?9 is not empty. Moreover, 
even though Ln$ might contain elements that are not functions (they are just Radon 
measure in such a case), its definition actually generalizes the usual logarithm. Indeed, 
for i9 G i?3-c + L 2 (0, T; Vb) we have (see [321 Rem. 4.3]) 

ln^GLn?? whenever $ > a.e. in Q and In $ G L 2 (Q). (4.4) 

Furthermore, the generalized logarithm is related to the theory of subdifferentials as 
follows. We define the function \1> : L 2 (0,T; V ) -> (— oo, +oo] by 

( v ) : = / rj,( v + 7?jc) for u G L 2 (0, T; V ) (4.5) 

being understood that the integral is infinite if ip(v + i?jf) ^ -^ X (Q)- Then, \1/ turns 
out to be convex proper and lower semicontinuous on L 2 (0, TjVb), so that its (possibly 
multivalued) subdifferential : L 2 (0,T;V ) ->■ L 2 (0,T;V *) is well-defined. Precisely, 
we have (see [121 Rem. 4.4] for details) 

Ln0 = 0tf(0-05c) for every d G ^ + L 2 (0, T; Vb). (4.6) 

At this point, we can state the generalized version of problem (I3.24p -f l3.26p as follows. 
We look for a quadruplet (tf,C,X,£) satisfying f l3TT5]) . fl337j) -f l3~TB|) . and 

C G -ff^O, T] Vq) and C G Lntf 

3s (C(*) + - «A0(t) = /(t) in Vb* for a.a. t G (0,T) 

<9 t X- AX + £ + c/(X) = A'(X)tf and £ G /3(X) a.e. in Q 
C(0) = ln^ and X(0) = X ■ 

The following result holds 

Theorem 4.2. ^sswme both (J33D-(J32D and (13T8D - d3~T3l) /iota and that k > 0. TTien, 
problem (14. 7p - (14.101) /ias a unique solution £, X, £) satisfying (13.151) ared (I3.17p - (l3.18p . 

Moreover, we have that 

lntf G L°°(0,T;if) and C = ln#- (4.11) 

Proof. We first prove uniqueness. Our proof closely follows [51 Sect. 5]. However, we 
repeat at least a part of the argument for the reader convenience. We first observe that, 
if CjX, £) is a generalized solution, by integrating (14. 8 p in time we obtain 

C + A(X) - kA(1 *■&) = ln^o + A(X ) + 1 * /. (4.12) 

Now, we pick two solutions Q, X$, i = 1, 2, and set for convenience 

0:=0i-# 2 , C:=Ci-C 2 , X:=Xi-X 2 , and £ := & - 6 . (4-13) 



(4.7) 
(4.8) 
(4.9) 
(4.10) 
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Now, we write ( 14.1 2\i for both solutions, take the difference, and test it by dQ t in the 
duality L 2 (0,T; V^*)-L 2 (0,T; V ), where t e (0, T) is arbitrary and Q t is the characteristic 
function of the interval (0, t). We observe that all the terms but one are in fact integrals. 
Namely, we have 



(C,^) + - / |V(l*tf)(i)| 2 = - / (A(X X ) - A(X 2 )) tf. (4.14) 



At the same time, we write ( 14. 9ft for both solutions, take the difference, multiply the 
resulting equality by X, and integrate over Q t . We obtain 

\ [ \x(t)\ 2 + [ ivx| 2 + / ex 

(A'(Xi) i?i - A'(X 2 ) tf 2 )X - / (^(Xi) - </(X 2 ))X. (4.15) 



Finally, we add ( I4.15P to (14. 14j) and proceed exactly as in [S]. However, let us point out 
that, in view of Taylor's expansion and the Holder and Sobolev inequalities, we have 



{-(A(X0 - A(X 2 ))tf + (A'(*i)#i - A'(X 2 )tf 2 )x} 
tfi{A(X 2 ) - X(Xx) - X'(X 1 )(X 2 - X 1 )} 

h 

+ / MX(Xi) - A(X 2 ) - A'(X 2 )(X 2 - Xi)} 
<c/ (^i + ^ 2 )|X| 2 < C / 11^(5) +i? 2 (a)|| 4 ||X(s)||4||X(a)|| 2 d3 



o 



<c / IliJiW+^WIIvOlVXWIIj+IIXWII^IIXWIIflds 



o 



^7/ |VX| 2 + c/ (\\Ms)\\ 2 v+\\Ms)\\ 2 y)\\X(s)\\ 2 H ds. 







Therefore, we derive that 



(C,^ t > + f / |v(i*tf)(t)| 2 + ^ / |x(t)| 2 + / |vx| 2 + / ex 

<c r(l + ll^( s )lly + ll^(^)Hv) Ms)\\ 2 Hds+ 1 - [ |VX| 2 . 



As the last term can be easily controlled by the left-hand side, what remains to observe 
is that all the terms on the left-hand side are non negative. The integral containing £ is 
non negative by monotonicity. Let us deal with the duality term. For % = 1,2 we have 
Q G Ln#;. Due to (14.31) . this means that 

(0, Vi - 0A + f V(^) < f ij>(vi) for every v i E$ i + L 2 (0, T; V ) and i = 1, 2. 
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Singular limit of a integrodifferential system 



Now, we choose the admissible functions V\ = #1 — $Q t and v 2 = $2 + $Ctj we sum U P 
and split the integrals. We have 



As v\ = §2 in Qt and v\ = $1 in Q\Q t and similarly for V2, all the integrals are finite and 
cancel out. We deduce that (£, — fiQt) < 0, i.e., what we wanted to prove. At this point, 
we can apply the Gronwall lemma (see, e.g., pp. 156-157]) and obtain, in particular, 
that X = and V(l * $) = a.e. in Q. As 1 * $ is Vo-valued, this implies that 1 * d = 
a.e. in Q, whence also d — a.e. in Q. All this means that $1 = ?9 2 and X\ = X^. By 
comparison in (14.121) and (14. 9p . we conclude that Ci = C2 an d ^1 = (2 as well. 

Once uniqueness of the generalized solution is proved, we can easily conclude. Indeed, 
our assumptions allow us to apply Corollary 13.21 Hence, a solution exists in the strong 
sense, i.e., satisfying the regularity requirements (I3.15p -( )3.18p . On the other hand, such 
a solution is also a generalized solution due to (14. 4p . Finally, it satisfies (14. lip . □ 



The argument we follow for our first proof uses compactness and monotonicity methods. 
So, we start estimating. However, we often proceed formally for the sake of simplicity. 
The correct procedure could be based on performing similar estimates on the solution of 
some approximating problem. One approximation is constructed in [5] and depends on 
the parameter e: the solution is smoother than the solution to the problem we are dealing 
with (actually the limit as e \ keeps such estimates). Furthermore, in order to simplify 
the notation, we often avoid the subscript r (on the solutions) during the calculation and 
restore it just at the end of each estimate. 

First a priori estimate. We would like testing (I3.19P by 



in the duality V^*-Vo and integrate over (0, t), where t £ (0,T) is arbitrary. However, 
we proceed formally, as just said. In particular, we behave as the logarithmic term were 
smoother. At the same time we multiply (I3.20p by d t X T and integrate over Q t (see (13.341) ). 
Finally, we sum up and remark that the terms containing d t X T partially cancel. Hence, 
by avoiding some subscripts in the notation for a while and adding the same integral 




5 Proofs of Theorems 13.31 and 13.4 



u T := i!} T 



(5.1) 
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Jo M^)| 2 to both sides for convenience, we obtain 

/ / dt(hx#)# x + [ (n \Vu\ 2 + (k T * Vu) • V 

l^l 2 + ^ / |VX(t)| 2 + / £(X(t)) + / |X(t)| 2 



/ dt\{X)ti K - f (K V# K +(k T *W x ))-Vu+ [ fu 

JQt JQt V J JQt 

+ I (|X(t)| 2 -a(X(t)))+ f ** + \ I |VXo| 2 + / (P + a)(X ). (5.2) 



2 . 

Note that we have •§ > and can make use of the chain rule for subdifferentials. Now, we 
recall that /3 is nonnegative (cf. (13.41) ) and treat each of the non-trivial terms, separately. 
We integrate the second integral on the left-hand side by parts with respect to time and get 



ft(lntf)^ = / # x (t) lntf(i) - / 7? M (0) In tfo - / hi W M 
< / Mt) - / Mt) hT&(t) + [ |ln#| \dtti % \+c. 



Hence, by recalling Proposition 14.11 and observing that r — d* ln + r > (r/2) — c for every 
r > 0, we deduce that 



•&{t)- / t (lntf)tf:K 
iQt 

> / U(t) -0*ln + 0(t) + In" i?(t)) - / |ln^||a t ^|-c 
in ^ ' JQt 

> [ (l#(t) + $Jn-l}(t)) - [ IIb^II^jcI-c. 
Jn yz J JQ t 

On the other hand, we notice that | lnr| < c((r/2) +0* In - r) for r > 0, so that 

I M \dt& K \ < c [ ((0/2) + A In" 0) 



and the last integral can be treated on the right-hand side via the Gronwall lemma since 
dt&x e L\0, T; L°°(fi)). Next, thanks to (IX27D and to (1X321 . we infer that 



/ (k |Vm| 2 + (fc T * Vu) ■ Vu) > f \Vu\ 2 >^[ 
JQt v y ig* ^ io 



i 

\\u(s)\\yds 



where C is the constant of ( I3.32p . Now, let us deal with the right-hand side. With the 
help of (13.271) . the Young theorem (13.371) . and ( 13.81) . we immediately have 

- / (k V^+ (k T * W k )) ■ Vu+ f fu<c + cf ||u(s)||^ds. 

JQt V J JQt io 

Moreover, we observe that 

dt\(X)0 x = f \'{x)d t xd n <- ! \d t X\ 2 + c[ \X\ 2 + c 

.In. 4 In. .In, 
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Singular limit of a integrodifferential system 



since |A'(r)| < c(l + |r|) by (13. 5p and < < by Proposition 14.11 Finally, we observe 
that (I3.5P also yields |cr(r)| < c(l + r 2 ) for every r and deduce that 

(\X(t)\ 2 - a(X(t))) < c + c [ \X{t)\ 2 <c + cf (|X | 2 + / 2X(s) d t X(s) ds 
! Jn Jn^ Jo 

<c+- \d t X\ 2 + c / \X\ 2 . 

4 J Qt JQt 

By combining all the estimates we have derived with (I5.2p . applying the Gronwall lemma, 
and owing to the Poincare inequality f!3.33j) . we conclude that 

||^T||L°°(0,T;L 1 (^))nL 2 (0,T ; y) + ||ln || L oo (0,T;L 1 (Sl)) + || X T || #1 (0,T;_H')nL°° (0,T; V) < C (5.3) 

besides an estimate for /3(X T ) in L°°(0,T; L 1 ^)). 

Second a priori estimate. We write (13.201) in the form of a nonlinear monotone elliptic 
equation, namely 

- AX T + £ T = X'(X T ) $ T - d t X T - a'(X T ) and £ r G (3{X T ) (5.4) 

and notice that each term on the right-hand side of (15.41) is bounded in L 2 (0,T; H) 
by (13. 5p and (15. 3p . Concerning the first term, notice that X T and d T are bounded in 
L°°(0,T;L 4 (fi)) and L 2 (0, T; L 4 (f2)), respectively, due to the Sobolev inequality ( 13T301) . 
Then, a quite standard argument (formally test (15.41) by either — AX T or £ r in order to 
estimate both of them and then use the regularity theory for elliptic equations) yields 

IMU 2 (0,T;W) + ||£t||l 2 (0,T;//) < C. (5.5) 

Third a priori estimate. We want to estimate d t $ T in L 1 (0, T; W~ 1,q (Q)) for some 
q > 1 satisfying L 2 (Q) C W -1 ' 9 ^), and the choice q = 4/3 will work. Therefore, by 
proceeding formally, we take any v G W 1,4 (Vl) satisfying \\v ||jyM(n) < 1 and test (I3.19P 
written at almost any time t by i9 T (t) v, which is a good test function belonging to Vo, due 
to (I3.30I) - (I3.31I) and (15.31) . We obtain (the first integral is intended as a duality pairing) 



dt-& T {t)v = {f(t)-\'(X T (t))d t X T (t)}# T (t)v 
Jn 

- [ «oW T (t) ■ V(0 T (t) V) - [ (k T * W T )0) ■ V(0 T (t) v). (5.6) 



We simplify the notation by dropping the time t (and the subscript r, as usual) for a 
while and estimate each term on the right-hand side of ( 15.61) . separately. We account for 
the Lipschitz continuity of A' (see (I3.5P ) and for the Holder and Sobolev inequality (13.301) . 
Moreover, we observe that ||f ||oo < c thanks to (I3.3ip . We have 



{/ - \\x)d t x}$v < \\f\Wh\\v\U + ||A'(x)|| 4 ||d t x|| 2 ||tf|| 4 |M 

<c||/|| 2 ||tf|| 2 + c(l + ||X|| 4 )||<9 t X|| 2 | 



'DO 



k W • V(dv) < k HWIIilMU + «o||Vt?|| 2 ||i?||4||Vi;||4 

n 

2 



< c \\vni+c\m\2\ 

f (k T *VA)-V(Av) < ||Av*V^|| 2 (||V7?||2||t;|| o + ||i?||4||Vt;||4) 

Jn 

<c||£; r *W|| 2 (||W|| 2 + ||tf||4). 
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We term C the maximum of the values of the above constant c's, for clarity. Then, we 
first collect the estimates just obtained and (15.61) . Finally, we take the supremum with 
respect to v G Wq (Q) under the constraint |M|wM(n) < 1. We conclude that 

ll^r(*)||w-M/3(n) < C<t> T {t) for a.a. t G (0,T) 
where T : (0, T) — > R is defined by 



Mt) ■= \\f(t)h\K(t)h + (i + \K(t)\U)\\d t x T (t)h\K(t)\U + \\vMt)\\l 

+ ||V7? T (t)|| 2 ||7? T (t)|| 4 + || (Av * W T )(t)|| 2 (||W r (t)|| 2 + ll^(t)|| 4 ). 
Therefore, the estimate 

II^T|Ui(0,T;W-M/3(n)) < c (5.7) 

follows once we prove that <ft T is bounded in L 1 (0,T). In view of the previous estimates 
(15.31) . (I5.5p . and of the Sobolev inequality, we see that the only trouble could come from 
the term containing the convolution. By owing to the Young theorem (see (I3.37P ) and to 
the Sobolev inequality once more, we have 

T ||(£v * W T )(t)|| 2 (||W T (t)|| 2 + ||0 T (t)|| 4 ) dt 

< c\\k T * W r 1 1 L*(P,T;H)\Wt\\l*(0,T;V) < 4 K ||x,i(0,T) \Wr || L*(0,T;V) ■ 

Now, we recall (13.271) 9 and (15.31) and conclude that r is bounded in L 1 (0,T). Therefore, 
(15.71) is established. 

Fourth a priori estimate. By testing (13.191) by any v G L 2 (0,T; Vq) and integrating 
over (0, T), we deduce that 

/ (d t \n$ T (t),v(t))dt < cM T \\v\\ L 2 {0>T . v) for every v G L 2 (0, T; V ) 
Jo 

where we have set 

M T := «o||V7? T ||£a( 0l T;fl) + ||AV * Vl9 r || L 2 (0iT;H) + \\f\\ L 2(0,T;H) 

+ (l + ||X r |U° o (0,T;£ 4 (n))) || 9 t X T \\ L 2 (0,T;H) ■ (5.8) 

Thus, the estimate 

||<9 t lntf T || i2(0)T .v-) < c (5.9) 

follows whenever we prove that M T < c. So, let us examine each term of (15.81) but the 
third one, of course, by accounting for (15.31) . The first and last ones are bounded by (I3.27P 
and (13.51) . For the second term, we use the Young inequality (13.371) and (I3.27P and obtain 

||AV * V^ r ||L2 (0iT;jH -) < \\k T \\ L i {0tT) \\V'& T \\L2(0,T;H) < C. 

Hence, (I5.9P is established. 
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Singular limit of a integrodifferential system 



Convergence and conclusion. From estimates ( 15. 3p and (15.91) we derive the following 



convergence 

$ T $ weakly in L 2 (0, T; V) (5.10) 

u T u weakly in L 2 (0, T; Vb) (5.11) 

X T ->X weakly star in H\0,T;H) nL°°(0,T;V) D L 2 (0,T;W) (5.12) 

£ r £ weakly in L 2 (0, T; if) (5.13) 

mtf T C weakly in H\0, T; V *) (5.14) 



for suitable functions X, £, (, and w := 'd—fin, possibly for a subsequence r = r n \ 0. We 
are going to show that the triplet ($, X, £) is a solution to problem (I3.24p - ()3.26p . Once this 
is proved, the whole family ($ T , X T , £ r ) converges (in the above topology) to ($, X, £), since 
the solution to problem ( 13.24l) -( j3.26p is unique. Now, we recall Theorem 14.21 and observe 
that it is sufficient to show that the quadruplet ($, (, X, £) is a generalized solution, i.e., it 
solves (I4.7p -( 14TT01) . So, we just prove this fact. The regularity requirements ( 13.151) . ( 13.171) - 
(I3.18p . and (14. 7p are already clear. First, we observe that (15.121) and (15.141) imply at least 
weak convergence in C°([0, T]; H) and C°([0, T];Vq), respectively, whence ( and X satisfy 
the Cauchy conditions ((0) = ln-#o and X(0) = X . Furthermore, as the embeddings 
W C V, V C L^(Q), and V C H are compact, we can apply, e.g., [T7J Sect. 8, Cor. 4] 
and deduce that (I5TT2]) . (IBTTUD and (J5ZZD imply 

X r -^X strongly in L 2 (0,T;V) n C°([0, T]; L 4 (fi)) (5.15) 
i? T i? strongly in L 2 (0, T; (5.16) 

Moreover, we can even assume that 

X T X and $ r -> $ a.e. in Q. (5.17) 

In particular, by noting that ( 13.51) implies Lipschitz continuity for A' and a' and the 
estimate (via Taylor's formula) 

|A(X r ) - A(X)| < c(l + |X|) \X T - X| + c\X T - X\ 2 

we infer that 

A'(X r ) ->■ A'(X) strongly in C°([0, T];L 4 (tt)) and a.e. in Q (5.18) 
a'(X T ) ct'(X) strongly in C°([0, T]; L 4 (fi)) and a.e. in Q (5.19) 
A(X r ) A(X) strongly in C°([0, T]; H) and a.e. in Q. (5.20) 

Next, we observe that (I5.10p implies that A-# T converges to Ad weakly in L 2 (0,T;V^*). 
By (I3.22p - (l3.23p and the Young theorem (see (I3.37P ). we infer that 

(k + 1 * fcr) * ^ k*M weakly in L 2 (0, T; V *) 

and we can take the limit in the integrated version of (I3.19p . Namely, we obtain 

C + A(X) - k * Atf = 1 * / + ln$o + A(X ) (5.21) 

and we conclude that the limit functions we have constructed satisfy (14.81) and 

d t X - AX + £ + o-'(X) = A'(X) a.e. in Q. 
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Thus, it just remains to prove that £ G /3(X) a.e. in Q and ( G Ln$. The first claim 
immediately follows by applying, e.g., [21 Lemma 1.3, p. 42] in the framework of maximal 
monotone operators in L 2 (Q) by accounting for (15.131) and (15.151) . On the contrary, the 
second claim needs much more work. We use the framework of the maximal monotone 
graphs in L 2 (0, T; Vo) x L 2 (0, T; Vq) and consider the subdifferential d^f of the function ^ 
given by (14.5)) . which is related to the multivalued operator Ln by (14. 6p . So, we have 
to prove that ( G d*&(u) by starting from ln$ r G d^(u T ) for r > (see (I4.4[) ) and the 
weak convergence given by (15. lip and (I5.14p . It is well known (see once more, e.g., [21 
Lemma 1.3, p. 42]) that a condition that allows to conclude is the following 

limsup (ln$ T , w r ) < ((,u), that is, limsup (ln$ r — (,$ T — #) < 0. (5.22) 

In order to prove (15.221) . we consider the integrated version of (13.201) and equation (I5.2ip . 
We take the difference and test the equality obtained by d T — d. We have 

(lntf T -CA-tf) = -ko / V(l * (t? t - 0)) • V(0 T - i?) 

- f (l*k T * W r - k' * W) • V(tf r -■&)-( (A(X T ) - A(X)) (0 T - i?) 

= -(/€ + «(,)/ V(l * (0 T - 7?)) " V(tf r - 7?) 



((l*fc r -«o)*W T ) -V(tf T -tf)- / (A(X T )-A(X))(0 T -0). 

Now, the first integral of the last chain is nonnegative and k + k' q = k > 0. Next, the 
last integral tends to zero by (15.161) and (15.201) . Finally, the middle term is estimated by 
the Holder inequality, the Young theorem, and (15.31) this way 



/ ((1 *Av -«(,)* V0 T ) .V(tf T -#) 
Jq 



< ||(1 * k T - K? Q ) * W r || L 2 ( o )T;if) ||V(tf r - &)\\v (0 ,T;H) 

< ||1 * k T - Ko||l,i(o,Dl|V^ T ||l ( a(o j T;H)||V(7? T - 0) \\v*(p,T;H) 

< c\\l * k T - k' \\ l i(o,t)- (5.23) 
Hence, by recalling (I3.22p . we conclude that (15.221) holds. 

Proof of Theorem 13.41 By arguing as in the last part of the previous proof, we consider 
the integrated version of the equations for temperature and test the difference by $ r — 
However, in the present situation we already know that ( = Iwd and can integrate over 
Qt rather than Q. Thus, a quite similar calculation yields 

/ (ln# T - ln0)(0 T -$) + k I V(l * (0 T - #)) • V(tf T - 0) 

((1 * k T - k' q ) * W r ) • V(0 T -■&)-! (A(X r ) - A(X)) (0 T - 0). (5.24) 

Qt JQt 



20 



Singular limit of a integrodifferential system 



At the same time, we multiply the difference between ( I3.20p and ( I3.25P by X T — X and 
integrate over Q t . We obtain 

\ \ \(X T -X)(t)\ 2 + I |V(X r -X)| 2 + f (£ r -£)(X T -X) 
z Jn JQt JQt 

(A'(X r ) tf T - A'(X) 0) {X T - X) - [ (a'(X T )-a'(X))(X T -X). (5.25) 



At this point, we sum (15.251) to (I5.24p . and it is clear that all the terms on the left-hand 
side are nonnegative. Thus, we estimate each term on the right-hand side. As far as the 
term containing the convolution kernel is concerned, we can repeat the argument that led 
to (15.231) . Thus, we have 

- / ({l*k T -K' )*VA)-V{$T-#)<c\\l*k T -K' \\ L l (0!T) . 
JQ 

The sum of all the terms involving A and A' can be first transformed and then estimated 
as follows (we use the Taylor formula and (13.51) . besides standard inequalities, as usual) 



tf r {A(X T ) - X(X) - \'(X T )(X T -X)}- / ${\(X) - \{X T ) + \'(X)(X T - X)} 

JQt 

<c I (tf r + tf)|X r -X| 2 
JQt 

<c [ \\# T {s) + ${s)\U\\X T {s)-X{s)\U\\X T {s)-X{s)\\ 2 ds 

JO 

< c [ \\A(s) + #(s)\\ v (\\X T (s) - X(s)\\ H + ||V(X T (s) - X(s))\\ H ) \\X T (s) - X(s)\\ H ds 
Jo 

\ j \V(X T -X)\ 2 + c f\l + \\Ms) + m\\ 2 v)\\Ms)-X(s)\\%ds. 

z JQt JO 



< 



As the integral involving a' can be treated in a trivial way due to (13. 5p . we can apply 
the Gronwall lemma (see, e.g., [SI p. 156]) and infer that the left-hand side of (13.281) is 
bounded by 



c||l*fc T -4||xi(o,T) exp(c j (1 + \\Mt)+$(t)\\v)dt). 
As the last integral is bounded by a constant thanks to (15. 3p . inequality (13.281) follows. 
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